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Abstract 

We study the models with the Majorana neutrino masses generated radiatively by two-loop 
diagrams due to the Yukawa and effective couplings along with a scalar triplet 

A, where p is a doubly charged singlet scalar, Ir the charged lepton and W the charged gauge 
boson. A generic feature in these types of models is that the neutrino mass spectrum has to 
be a normal hierarchy. Furthermore, by using the neutrino oscillation data and comparing with 
the global fitting result in the literature, we find a unique neutrino mass matrix and predict the 
Dirac and two Majorana CP phases to be l.dOvr, l.llvr and IATtt, respectively. We also discuss 
the model parameters constrained by the lepton flavor violating processes and electroweak oblique 
parameters. In addition, we show that the rate of the neutrinoless double beta decay (On/3/3) can 
be as large as the current experimental bound as it is dominated by the short-range contribution 
at tree level, whereas the traditional long-range one is negligible. 
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I. INTRODUCTION 


Although the data 
carry nonzero masses 


rom neutrino experiments have implied that at least two neutrinos 




the origin of these masses is still a mystery. Apart from the 


mass generation of Dirac neutrinos given by the Yukawa couplings with the existence of 


right-handed neutrinos (zz/j), seesaw mechanisms with type-I [6|-[l0(], type-II 
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and 


type-III 18| can generate masses for Majorana neutrinos by realizing the Weinberg operator 
Ll) at tree-level, where $ and Ll are the doublets of Higgs and left-handed lepton 
fields, respectively. In these scenarios, either heavy degrees of freedom or tiny coupling 
constants are required in order to conceive the small neutrino masses. On the other hand. 


models with the Majorana neutrino masses generated at one-loop 


and higher loop 


25 


Q,Q, 


two-loop 


2l|-l24| 


28| diagrams have also been proposed without introducing z/^. Due 


to the loop suppression factors, the strong bounds on the coupling constants and heavy 
states are relaxed, resulting in a somewhat natural explanation for the smallness of neutrino 
masses. 


Among t 
models |23|, 


le loop-level mass generation mechanisms, there is a special type of the neutrino 
24| in which a doubly charged singlet scalar p : (1,4) and a triplet A : (3,2)^ 


under SU(2)i x U(l)y are introduced to yield the new Yukawa coupling pi^in with the 


mixing 


charged lepton in as well as the effective gauge coupling p^^W^W^ due to the 
between and leading to the neutrino masses through two-loop diagrams 23|. As 

this model is the simplest way to realize the pWW coupling, we name it as the minimal 
two-loop-neutrino model (MTM) 23|. It is interesting to note that p^'^W^W^ can also be 
induced from non-renomalizable high-order operators 29l-l3l|. Although MTM can depict 


neutrino masses at two-loop level, the assumption on the absent of the L'^LA term makes this 
model unnatural. To solve this problem, one can simply extend MTM by adding an extra 
doublet scalar, which together with A carries an odd charge under an Z 2 symmetry 321. 
We call this model as the doublet two-loop-neutrino model (DTM). On the other hand, 
p^^W^W^ could be granted by inner-loop diagrams, such as those |27|, |28| with three-loop 
contributions to neutrino masses, in which the neutral particle in the inner-loops could be 
a candidate for the stable dark matter. 


^ The convention for the electroweak quantum numbers {I,Y) with Q = / -|- y/2 is used throughout this 
paper. 
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In this study, we will demonstrate that the neutrino mass matrix can be determined 
in these models by the experimental data. In particular, the neutrino mass spectrum is 
found to be a normal hierarchy. In addition, the neutrinoless double beta decay (Oz//3/3) 


is dominatec 


W the short-ran ge c ontribution at tree level due to the effective coupling of 


23l. I 24 J. I28l430l. I 33 I. l34j |. instead of the traditional long-range one. However, the 


neutrino masses in this type of the models are usually over suppressed as there is not only a 
two-loop suppression factor, but also a small ratio mi/v with the charged lepton mass mi and 
vacuum expectation value (VEV) v = 246 GeV of the Higgs held. Furthermore, the lepton 
havor violation (LFV) processes could also limit the new Yukawa couplings. To have a large 
enough neutrino mass, the mixing angle or mass splitting between the two doubly-charged 
states should be large, which inevitably leads to a signihcant contribution to the electroweak 
oblique parameters, especially the T parameter. We will calculate the neutrino masses in 
details and check whether there is a tension between these masses and the constraint from 
the oblique parameter T. 

This paper is organized as follows. In Sec H, we study the neutrino masses in the two- 
loop neutrino models. In Sec HI, the constraints on the model parameters from lepton 
flavor violating processes and electroweak oblique parameters are studied. We present the 
conclusions in Sec. IV. 

II. TWO-LOOP NEUTRINO MASSES 


In MTM, we introduce the scalars p : (1,4) and A = (A++, A+, A°) : (3,2) under 
SU{2)l X U{1)y- The relevant terms in the Lagrangian are given by 


-£ = -/i|(<I)t$) + M2(AtA) + A$(<I)t$)2 + A3(AtA)i(<l>t<I))i + A4(AtA)3(<I)t<I>)3 
Yah{Ll)aA{LL), + ^p(4),(£^), - /iA(4>t)2 + + Xp*A^^ + H.c. 


( 1 ) 


where <h = (d*’'', with <l>° = ( 4 ) 7 ^ -|- is the SM doublet scalar, the indices of a 

and b represent e, p and r, and the subscripts of 1 and 3 in the quartic terms stand for the 
SU(2) singlet and triplet scalars inside the parentheses, respectively. After the spontaneous 
symmetry breaking, 4) acquires a VEV of = \/2 (4*°), while the neutral component of A 
also receives a VEV va/\/2, generated via the p term. Note that by the global htting result 
of po = 1.0000±0.0009 35|, ua is constrained to be < 5 GeV, so that ~ 246 GeV is a good 
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FIG. 1; Diagrams for the neutrino mass generation, where the charged states can be replaced 
by Sf or when DTM is discussed. 

approximation. The k term in Eq. ([1]) can produce a mixing term between and 
resulting in two mass eigenstates Pi ,2 with masses Mi, 2 , respectively. We will set Yab and A 
to be zero since they have the tree-level and logarithmic divergent two-loop contributions to 


neutrino masses, respective’ 
introducing a new doublet 


y. These two coupling can also be forbidden in a natural way by 


32l | or a singlet scalar 33(] with an Z 2 symmetry or by replacing 


A by a higher multiplet, such as ^ : (5,2) without the discrete symmetry 3^. The scalar 
mass spectra of MDM are shown in Appendix Rll. 

We now calculate the neutrino masses from the two-loop diagrams of Fig. [T]in the t’Hooft- 
Feynman gauge. The neutrino mass matrix can be written as 

1 




— [^(a) + ^(b) + A(c) -b , 


( 2 ) 


(Ibvr^)^ V" 

where the integration results A(a), A(^b), A(c) and corresponding to the sub-hgures 

(a), (b), (c), and (d) in Fig. [H are given in Eqs. (B9)-(B11) in Appendix iBl respectively. 
Explicitly, we find that the contribution related to A(^a) dominates over the other three 
components. We note that if Mp is much smaller than Ma and the mixing angle 9 between 
them is small, this model approximately coincides with the effective theory involving the 
dimension-7 operator p(P^<I))(P^<h)<h<h discussed in Ref. [30]. 

DTM can be viewed as the extension of MTM by introducing a new doublet y = (y*", 
with y° = (yb -\- iyj) / \/2. This new doublet along with A carries an odd charge under the 


Z 2 symmetry 


32l |. This discrete symmetry can forbid the tree-level coupling L^LA to make 
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the two-loop neutrino mass generation more natural. The relevant part of the Lagrangian 
is given by 


- £ = -/i|(<ht<h) - + A$(<ht$)2 ^ + A4(<htx)(x^$) + Ml{A^A) 

+ [^PiWR)b - pA^^X^ + ^P*A^ + \p*A^x + + h-c." . 


(3) 


Since y does not couple to the SM fermions due to the Z 2 symmetry, the model is similar to 
the Type-I two-Higgs doublet model 37|. The doublet y can also have a VEV v^/\/2 = (x°) 


due to the negative mass term of x- We can dehne the mixing angle sin 7 = n$/ v'^ 

to characterize the scalar mixing matrix when only scalar doublets are taken into account, 
where v‘^ = v% + v‘^ + 2v\ with v = 246 GeV. Note that the A 5 term in Eq. (j3]) breaks the 
lepton number symmetry explicitly so that the dangerous Majaron can be avoided, while 
sizable values of A$ and A^ are needed in order to give the CP even neutral scalar masses 
and preserve the stability of the scalar potential. 

The VEV of A in this model is induced by the p term in Eq. ([3]), which is proportional 
to with hxed values of Ma and p. In the model, we have two singly-charged physical 
states Sf and S 2 besides the unphysical Goldstone boson originated from the mixing 
among and A^. The values of the mixing elements between the doublets and 

A are also proportional to ua like MTM. Moreover, the term \p*A^x ci-nd its hermitian 
conjugate provide another source for the — A^^ mixing apart from the k term, with 
the contribution to 9 approximately proportional to sin 27 . The results on the scalar mass 
spectra are given in Appendix El 2. 

The mechanism for the neutrino mass generation in DTM is similar to that in MTM. 
But, the main coupling related to is from the effective dimension-5 effective operator 
p(x*h)^. The formula for the neutrino mass matrix is given by 




2Cuimem£i 


v‘^c^ 


+^(6) + ^(c) + A(d) 


S2e 


K 


Xv 


P — A{al) + 77^(a2) + 


(a3) 


(4) 


where A(^ai) and A(j) with z = 1, 2 and 3 and j = b,c and d are dehned in Eqs. (B12)-(B15), 
respectively. In Eq. (01), there is a new contribution proportional to A, which is of 0{v/\/v). 
Note that the elements of the neutrino mass matrix in MTM are of 0{v\/v‘^). 

It is crucial that the above types of the two-loop neutrino mass generation, in which ^r^rP 
is the only source of the LEV, can lead to an interesting structure for the neutrino mass 
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matrix. The relative sizes among the matrix elements are determined by the combination 
factors of Cuimim(,i. Assuming that each value of Cu/ is at the same order, there exist 
interesting hierarchies for the mass matrix elements, given by 




(5) 


In particular, (Mj,)ee is much less than due to ml/rri^ ~ (1/3500^). In Refs. |38l-l43|. 

it has been shown that only the normal hierarchy for the neutrino mass spectrum can have 
the matrix textures in which {M,^)ee together with another matrix element is zero. Clearly, 
as the mass hierarchies in Eq. (|5]) naturally realize {M^)ee — — 0, both MTM and 

DTM predict the normal neutrino mass hierarchy 


Recall that in the standard parametrization 35|, l4^, the neutrino mixing matrix Vrmns 
is given by 


hpMNS — 


C12C13 

’S12C13 

SI3& 


1 0 

0 ' 

-S12C23 - Ci 2 S 23 'Sl 3 e*'^ 

C12C23 - Si 2 S 23 Sl 3 e*‘^ 

S23CI3 


0 

0 

i S12S23 - Ci 2 C 23 Si 3 e*'^ 

-C12S23 - Si 2 C 23 Si 3 e*'^ 

C23CI3 J 


0 

0 

gia 3 l /2 


( 6 ) 


where Sij{cij)=sin9ij {cos 6 ij) with 6 ij being the mixing angles, 6 is the Dirac phase, and 
021 and Q! 3 i are the two two Majorana phases. Note that one of the Majorana phases can 
be absorbed by the chiral fermion helds if there exists one massless neutrino. For given 
values of mass square splittings and mixing angles, there are only two solutions for the three 
CP phases of S, 021 and 031 , along with the lightest neutrino mass mo, to satisfy the mass 
hierarchies in Eq. ([5]). In particular, by using the central values of the global htting result 
for the normal hierarchy mass spectrum, given by 351 


2 012 = 0.308 ± 0.017 , sin^ 623 = 0.437l°;°| , sin^ ^13 = 0.0234l°;°°20 , 


sm 6^12 


Am2i 


we hnd that 


—u.uzci ’ (j.uuiy ’ 

= (jMtoiij X 10"^ eV , Aml^ = (2.43 ± 0.06) x 10“^ eV, 


(7) 

( 8 ) 


(i) : mo = 5.14 X 10 ^ eV , 5 = 0.60 vr , 021 = 0.11 tt , 031 = 0.53 tt , 

(ii) : mo = 5.14 x 10“^ eV , 6 = 1.40 tt , 0^21 = 1.11 tt , 031 = 1.47 tt . 


(9) 

( 10 ) 


Note that both solutions in Eqs. ([9]) and (fTOj) have the same value for tuq but different CP 
phases. It is interesting to see that the predicted Dirac phase 6 = 1.407r in (ii) of Eq. flTOj) 
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agrees well with that given by the global fitting result in Ref. 3^. Taking (ii) in Eq. flTOj) 
as the input parameters, the neutrino mass matrix is then given by 


M„ = 


( 


\ 


0 0 1.04e*i-937r 

0 2.42 2.32 

l_04e*i-93^ 2.32e*°-^°^ 2.79 




/ 


( 11 ) 


in unit of 10 GeV. Note that the empty values for [My)ee and (Mj,)e^ can be placed by 
some small non-zero values when any of the parameters in (ii) is under slightly shifting. 


III. CONSTRAINTS FROM LEPTON FLAVOR VIOLATION PROCESSES AND 
ELECTROWEAK OBLIQUE PARAMETERS 

In both MTM and DTM, as the coupling matrix elements Cab are the only sources of 
the LEV, the processes of £ ^ £' 7 ) with the tree-level (one-loop) contributions 

involving could give signihcant constraints on Cab- However, those on Cee and Ce^ can 
be ignored since they do not affect the tiny matrix elements (M^)ee and when we 

discuss the neutrino mass spectrum. Among the current experimental bounds, —)■ 

c'"7) < 5.7 X 10“^^ ^ is the most stringent one to limit of Cab- In particular, we can 
obtain 


la 


|2 / C 


M2 ^ M? 


< 


0.336\2 


TeV 


( 12 ) 


‘1 ^>'-'2 

It is obvious that the largest allowed value of |Ger|max from Eq. (IT^ depends only on Mi since 
se is of order 10“2. To account for the current experimental data on the neutrino masses as 
obtained in Eq. flTTl) . the matrix element (Mj,)eT should be around 1.04 x 10“^^GeV. As a 
result, we can use this value to check whether the mechanism of the neutrino mass generation 
can work, as shown in Fig. [2l The value of k is taken to be k < max(Mi, M 2 ), constrained 
by the perturbativity In general, a larger allowed value of k is more possible to give a 
correct value of {My)f.T- To obtain the right values for the neutrino masses, at least one of 
Ml and M 2 should roughly larger than 2.5 TeV. In Fig. [2^, {My)eT behaves approximately 
as an increasing function of M 2 due to the weak bound on C^' from the LEV processes, 
while in Fig. | 2]3 it is linearly proportional to M 2 . 
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(a) Mi<M 2 (b) Mi>M 2 



FIG. 2: Plots for with (a) Mi < M 2 and (b) Mi > M 2 . 


The experimental constraints from the /i —)■ e conversion could also give some hints on 


Ml. To illustrate the result, we pick out some of the experimental bounds, given by B 


Au 


7x 10 


-13 


47j, <7 X 10-1^ |48|, Bj^^ < 4.3 x < 45 ^ iq-u 


< 


For MTM and DTM, the dominant contributions come from 7 and Z penguin diagrams, 
which lead to 




2GFm\ 


pcapt 

^ A 


eALD(A) + Sa'^A) + sS' 


(n) ■ 


(13) 


where is the muon capture rate for the nucleus A, the coefficients Al and cor¬ 
respond to the dipole and vector contributions, and D{A) and VP’^(A) are the overlapping 
functions between e and /r (see Ref. for the details), respectively. Explicitly, we have 


At, = 


3\/2 


1927r2M2GF 


CiilC*i , gRViq) 


’IV 


SM, 


w 


27r2 9M2 




(14) 


i p i 

where Qq is the electric charge of the quark q and gRv{q) is the vector coupling with the 
quark q, mainly from the 7 penguin diagram as the Z diagram is suppressed by the charged 
lepton masses. Based on the valence quark model, one has the relations between g^^y and 
9rI, given by g^^l = 2gRy^u) + 9RV{d) and g^^^ = 2gjiyi^d)+9RV{u) |^- By taking Mp = ITeV, 
Cer = 0.33, and C^r = 0.0033, we hnd 


= 1.4 xl0-'^RS^, = 8.3x10 


-15 


= 9.8 X 10 


-15 


l-i^e 

?A1 


= 1.2x10-^^ 


, B;^^ = 7.5 X 10 


-15 


(15) 


which satisfy all the corresponding experimental limits. The improvement on the sensitivity 























of the fi — e conversion in the future will either detect the signal or put some more 

stringent constraint on the models. 

It is interesting to note that the neutrinoless double beta decay in our models can have a 
signihcant different feature from other models with radiative neutrino mass generations. In 
MTM and DTM, the short-range contributions to the decay dominate over the traditional 
long-range ones [23|, [2^, with the decay amplitudes proportional to {M^)ee- It is clear 


that the long-range parts can be safely neglected due to the small electron mass in (Mj,)ee, 
whereas the short-range ones are proportional only to the Yukawa coupling Cee- As a result, 
by calculating OiyjSlS, the upper limit on \Cee\ could be derived, despite of the fact that it is 
ignored when discussing the neutrino mass matrix. The half life for Oz//?/? is given by 54 155| 


rpOi/fSp _ \^LLL\2\ A/< |2a —1 

^ 1/2 “ 


1-Ma 


which leads to 


29| 


,LLL 


= mp {2C;^S2e[ 


V/\ Ml — Ml 


(16) 


(17) 


\/2 Ml Ml ’ 

where Gqi and \A4.^\ are the phase space factor and the matrix element for the hadronic 
sector, respectively, and is the coefficient, which is effectively related to the dimension- 
9 operator {uL'jpdL){uL'y^dL){eRe']^), dehned in Refs. [54, l55|, and nip is the proton mass. 
Note that the coefficient in Eq. (ITTjl has no explicit dependence on the electron mass. If 
one takes Mi = 1 and M 2 = 1.5 TeV in MTM, resulting in the maximal value of mixing 
I sin2^1 = 0.04, the upper bounds on \Cee\ for different target nuclei can be estimated as 


< 0.088 , < 0.067 , < 0.36 , 

< 0.096 , < 0.36 , < 0.13 , 


(18) 


by comparing with experimental upper limits 


56 


62| . When including the effect of A 


m 


DTM, a larger contribution to sin26* could make these upper bounds on \Cee\ to be around 
20 % lower. 

Combing the typical value of {M^)eT and the constraints from the LEV processes, at least 
one of Ml and M 2 should be heavier than around 2.7 TeV in MTM. On the other hand, to 
get Ml = (9(10^) GeV, M 2 needs to be much larger, at least 4 TeV. However, it is more 
difficult to get the value of M 2 less than 1 TeV with a large Mi, which means that to get 
Mp = C>(100)GeV, at least Ma > 4 TeV is required. Gonsequently, it is possible to detect the 
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FIG. 3: Allowed regions in Mi — M 2 plane for (a) MTM and (b) DTM, where k = 5max(Mi, M 2 ) 
and Cer = (C'eT-)max are used in MTM and DTM, while = 0.4, A = 2, A 4 = 1 and A 5 = 1.5 are 
taken in in DTM. Gray areas represent regions without enough neutrino masses, and blue regions 
located at Mi ~ M 2 are disallowed due to the mixing term between and 


signals of mainly throngh the pair prodnction of and 

with the same sign charged leptons in the hnal states at the LHC 


)he 

su 

30, 

63 


related resnlts in Fig. [3^. Similar conclusions have been also shown in Fig. 2 and 8.8 of 
Refs. 33, respectively, but they allowed some of the region with Mi M 2 ~ 1.5 TeV, 
which is forbidden in this paper. In DTM, the neutrino masses could be lifted up more 
easily by using a sizable A as well as sin0. Moreover, there is a new contribution to the 
neutrino masses from the dimension-5 operator p^^x^ DTM instead of the dimension-7 
one p[(D^<l>)<l>]^ in MTM. As an example, if we take = 0.4, A = 2, A 4 = 1, and A 5 = 1.5, 
along with the same values of k and Cer in MTM, one hnds that M 2 < 400 GeV (> 2 TeV) 
for Ml > 550 GeV (< 400 GeV). The relevant result is displayed in Fig. [3b. 

Finally, we briefly discuss the effects of the electroweak oblique parameters S and T in 
our models. First of all, in MTM and DTM, we hnd that the typical value of S is of order 
10“^, which is lower than the current experimental sensitivity. For the T parameter, the 
mixing between A and <1> gives a logarithmic divergent T due to the non-unity of p, but 
this part could be ignored when v 5GeV. In this case, the main contribution to T is 
given by the mixing between A and p, denoted as T(p_A). It is basically a negative value 
whose magnitude is limited due to the small mixing angle |s6i| < 0.02. For example, taking 
Ml = 2 TeV, M 2 = 4 TeV, = 0.5GeV, and k = 5 M 2 , it only leads to T(^p-A) = —5 x 10“^. 
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However, in DTM the mixing angle |s 6 )| can be enhanced by A, which is independent of 
va- Using the input values for the above parameters, and A = 3, we get T(^p-A) = —0.001. 
Meanwhile, the mixing between $ and y can also provide a sizable value to the corresponding 
parameter For example, T(^<s>-x) = 0.04 with Sa = = 0.4, A 4 = 1, and A 5 = 1.5, 

which still fulhlls the experimental bound —0.02 < AT < 0.12 at 1.5 — 1.8 a conhdent 
level 35(] . The relevant formulae for T are summarized in Appendix O 


IV. CONCLUSIONS 


We have studied the two models of MTM and DTM, in which Majorana neutrino masses 
are generated radiatively by two-loop diagrams due to the Yukawa pip^iR and effective 
couplings. We have shown that the lepton violating processes, in particular, 
^ e+ 7 can give stringent constraints on the new Yukawa coupling Cui. By combing 
with the perturbativity condition for the coupling k, the light neutrino mass element {M^)eT 
can limit the allowed values of Mi and M 2 . In particular, we have found that only Mi and 
M 2 with TeV scale can lead to the correct sizes of the neutrino masses. We have illustrated 
that the normal neutrino mass hierarchy is a generic feature in these two-loop neutrino mass 


generation models. Moreover, by using the central values of t 


le neutrino oscillation data and 


comparing with the global htting result in the literature 


35j, we have obtained the unique 


neutrino mass matrix in Eq. fllll) and predicted the Dirac and two Majorana CP phases to 
be l.dOvr, l.llvr and IATtt, respectively. Finally, we emphasize that the neutrinoless double 
beta decays can be very large as they are dominated by the short-distance contributions at 
tree-level, which can be tested in the future experiments and used to constrain the element 
of 
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Appendix A: Scalar mass spectra in MTM and DTM 


MTM 


The non-self-Hermitian terms in Eq. ([T]) can be expanded as fellows: 

v 2 

p*A^ = 2p* A++A°-^(A+)^ 


(Al) 

(A2) 


After obtaining the explicit forms of the scalar potential, we can write down its tadpole 
conditions 


—p\vq, + A<i,n| — V2pv^VA = 0 , 
MinA-/i^ = 0, 


which give 


Va = 


h vl 


V 2 MI ■ 

The mixing matrices of CP odd neutral and singly charged states are written as 

+'2 _+/ 

h/q bo 



M 


A ) 


(A3) 

(A4) 

(AS) 

(A6) 


where = 2vA/vq, and = \/2vA/vq>. The masses of the neutral CP odd state A and 
singly charged states are given by = Ma/c^ and Ms = Ma/c|^ respectively. 

The doubly-charge mixing matrix is given by 


Ml 


where 


A/g± = 


KUa 


One can easily diagonalize Eq. flA7p . leading to the mixing angle 6 


t2e — 

and eigenvalues of the eigenstates Pi ,2 

Ml = Ml - 


2Mlf 

Ml-Ml ’ 


’0 


(A7) 


(AS) 


(A9) 


Ml = Ml + 




2 (Mi - Mi) . 


1 - 2s 


(AlO) 
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2. DTM 


The related operators with x in Eq. ([3]) can be written as 

A+ 

\/2 

r A^ 

p*A^X = P~~ A++<h°x°-+ A°$+x’' 

L v 2 

The minimization conditions are given by 

-/r|n<i, + A$n| - ^(A 4 + A 5 )n$n^ - -^v/^v^ = 0 , 

-t4'’x + - ^(-'i + '\5)t>»4 - = 0 > 

MIva - = 0 ■ 


(All) 


(A12) 

(A13) 

(AM) 


It is convenient to define v = ynl + and = v^/v. The singly charged and CP odd 
neutral mass matrices are both 3x3. In the diagonalization, we use the relation ^ s'j^. 
The transformation matrices, Vj and V±, of CP odd neutral and singly charged states can 
be presented by a set of small quantities Cjj and eC, given as 


where 



/ 




/ 


, \ 



- Sry 

ei3 



— s^i 

ei3 

C/ = 



£23 

,^± = 



^3 


\ *"31 

^32 

1 J 


\ 

J 

^32 

1 / 


(A15) 


and 


^32 — 


Ml 


C2 a 


Ml - S2^ ^ ’ 


^31 — tp , 

_ ^ ,^27 Ml 


^23 — i.—S'r 


C27 Ml 
2s^Ml-X,v^’^ ^ 


(A16) 

(A17) 

(A18) 

(A19) 


/ ^ Ml C27 ■ 

Mi-i(A4 + A5)fi2s27 

^31 ~ ^0 ) 

C27 __ 

2c^ Ml — ^(A4 + As)"!)^ 
C27_iW|_ 

2s^ Ml — ^(A4 + As)!)^ 



t'. 


'0 ’ 



(A20) 

(A21) 

(A22) 

(A23) 
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The mass eigenvalues for two CP odd states are and M\^ = while those 

for singly charged states = A(A 4 + and = M\. Finally, the mixing angle 

between the doubly-charged states is given by 


tie — 


\/2kva + Ac^s^n| 


Mj-Ml 


(A24) 


For the mixing of CP even neutral states, one can focus on the 2x2 mixing matrix between 
$ and X) given by 


I 2A$n| -(A 4 A 5 )n$nx 

y -(A4 

Rotating the states = Cah — SaH and xr = -|- CaH, we obtain two mass eigenvalues 

Mfi and Mr to be 



— ( 2 A$ — (A 4 -|- X^)t^ta)v\c!p(^ , (A26) 

Mr = (2A$-|-(A 4 -|-A 5 )^)n|c||c^ . (A27) 


If one takes m/j = 125.7GeV [3^ as the state found at the LHC [ 6 a, Mr can be obtained 
as a function of A 4 -|- A 5 and a. 


Appendix B: Two-Loop Neutrino mass 

1. MTM 


We calculate the neutrino masses by using the t’Hooft-Feynman Gauge. The relevant 
vertices corresponding to the top vertex in Fig. UK are —pA and — ■|p++*(A+)^, 

which yield 


— 


ivA fCu'memp's^ 


(167r^)2 V v'^Cp, 


/kuaX m2 


111 


2/121 + h 


221 


+ sl{I_ 


112 


2/212 + / 


222, 


(Bl) 


with 

rl r^—yi j•l—x2 2 

Iijn= / dyi / dyi / dxi / dxi- - \og{yi[xiMl + X2M’^] + y2Xi{l - Xi)Mj‘^), 

Jo Jo Jo Jo t Xi 

(B2) 
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where we can define M[ = Mw and = Ms- For Fig. [Dd, we get 




i / y/2Cuim£m£i \ V^dw '^26 ^2/3' . ,^2 


(167r2)2 V / 2 2 2 


(M? - , 


(B3) 


with 


Ab) ^ 


r\ f\ ri--2-^3 _2y,{2-x,-X2) 

/ dy2 / dyi / dx^ / 0x2 / axi- 

'0 -/o -/o -/o -/o 


(xi + a;2)(l - xi- X2)m. 


2 ’ 


(B4) 


where 


nil = yi{xiMl + X 2 MI + x^Ml,) + |/ 2 (a:i + X2){1 - Xi - X2)M' ^ . (B5) 

The amplitnde shonld be eqnal to ■ For Fig. [T]i, we have 

- -Mf = (I6~7ry s,,{m,m,,Ca.)in^ (J^i) - /<‘‘I(M4 )|(p£>.0 , (B6 ) 

rl r^-yi rl rl-x 2 _a 

= / dy 2 dyi / cix 2 / , (B7) 

where 


fnfd) = yiixiMi + a;2M^) + 2/2a;i(l - Xi)M^ 


In snmmary, A(a), 24(b), A{c)i and A(^d:) in Eq. (|2]) are listed as follows: 

A/f 2 r o '2 

7 l(a) Sf^fK- 


Ml 


Mf - Ml 


(-fill — -^112) H—^(-^121 


V 2 


J122 + 7211 


7212 ) 


„/4 

„/4 

C/3 


(7221 — I 222 ) + si/K cl{I 


111 


'121 — -'211 


+ 7; 


221J 


+sl{I 


112 


'122 ~ -<212 


I212 + 7 


222 ) 


A{b) — A(c) — 2s^,kM^(7^^ 7^) , 


A(rf) = 2s%k 


j(d) 


Mf - Ml 


(B 8 ) 


(B9) 

(BIO) 

(Bll) 
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2. DTM 


If we expand the amplitudes A(ai), ^(a 2 ), and ^(as) up to 0{yi^/v), then the results are 
given by 


^(al) — 0 , ^(a2) — 0 , 


(B12) 


^(a3) ^"f^d 


2c 


'7^31/111 + 2c^st^e'^i{Ii2i + /211) + 2c^e^]^3(/i3i + /311) — 26325^/221 


+25^633(1231 + /321) + 26^57633/111 — 26757(63367 + e32 57)(/l21 + /2I1) 

+26757633 (/l31 + /311) + 2675^632/221 ~ 26757633 (/231 + /321) 


-+ 6 ^ 


26^57/333 — 26757(6^ — 5^)(/321 + /2I1) — 2675^/221 


A _ A _ "^1 "^2 1^2 / J t(P) „ / rib) . I T(.b)\ 

^{b) — ^(c) — 2S2e ^^m+31-'iv '^7^32-'5i + ^13-'52 I ’ 


V 


/!(,) = 2 s; 3 '—M^/W . 

V 


(B13) 

(B14) 

(B15) 


where we have used the same notations as those in Appendix [B]l except = Ms^ and 


M' = Ms,. 


Appendix C: T parameters in MTM and DTM 

The T parameter due to the mixing between p and A in MTM is given by 


T(p-A) = 




slF{MlMl^) + clF{MlMl^) - F{MI,MI^) 


2cisiF{Mt,M^ 


(Cl) 


The above result in Eq. fICip can also be applied to DTM with replacing by S^. In 
DTM, since the contribution from the mixing between $ and x should also be considered, 
we hnd 




167r5^M^ 


^a—'y 


M^) - G{Ml Ml)) + sl_JG{Ml, M^) - G(iV/^, M|)) 


+sl_^{F{Ml Ml) - F{Ml Ml)) + cl_lF{Ml, Ml) - F{Ml, Ml)) 
+FiMl,Ml: 


(C2) 
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where the functions F, K, and G are dehned by 


G{x,y) = F{x,y) + 4:yK{x,y) , 


(C3) 

(C4) 



X log X — y log y 


x-y 
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